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ABSTRACT 

We present XMM-Newton observations of NGC 891, a nearby edge-on spiral galaxy. 
We analyse the extent of the diffuse emission emitted from the disk of the galaxy, and 
find that it has a single temperature profile with best fitting temperature of 0.26 keV, 
though the fit of a dual-temperature plasma with temperatures of 0.08 and 0.30 keV 
is also an acceptable fit. There is a considerable amount of diffuse X-ray emission 
protruding from the disk in the NW direction out to approximately 6 kpc. We analyse 
the point source population using a Chandra observation, using a maximum likelihood 
method to find that the slope of the cumulative luminosity function of point sources 
in the galaxy is — 0.771° i 3 - Using a sample of other local galaxies, we compare the 
X-ray and infrared properties of NGC 891 with those of 'normal' and starburst spiral 
galaxies, and conclude that NGC 891 is most likely a starburst galaxy in a quiescent 
state. We establish that the diffuse X-ray luminosity of spirals scales with the far infra- 
red luminosity as Lx oc L a F1 7 j^ ' a7 , except for extreme starbursts, and NGC 891 does 
not fall in the latter category. We study the supernova SN1986J in both XMM-Newton 
and Chandra observations, and find that the X-ray luminosity has been declining with 
time more steeply than expected (Lx oc t~ 3 ). 

Key words: galaxies:individual: NGC 891 - X-rays: diffuse emission - galaxies: star- 
burst - supernovae:individual: SN1986J - X-rays: galaxies 



1 INTRODUCTION 

The environment of a spiral galaxy can influence its evolu- 
tion in many ways. Close interactions with nearby galaxies 
have been shown to trigger inspire star formation, and thus 
alter the interstellar medium (ISM) of the galaxy. Prop- 
erties of the ISM could also be affected by direct infall 
of material from the intergalactic medium (IGM), or, if 
the galaxy is in a group or cluster, by the ram pressure 
stripping by the IGM. However, in certain, outflows into 
the IGM could be a significant mode of the interaction of 
a galaxy with its surroundings. Bursts of star formation 
in the galaxy can, if sufficiently strong, lead to the out- 
flow of matter, which w ill be predominantly visible at Ha 
and X-ray wavelengths iShopb ell fc Bland-Hawthornlll998t 
iLehnert. Heckman fc Weaven^99l) . These outflows can be 
strong enough to escape the visible extent of the galaxy, but 
still be retained by the halo of the galaxy. In some of these 
cases, the burst of star formation would be strong enough 
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for the outflow to be able to could escape the galaxy and its 
halo, and appear as a superwind. 

A detailed study of the region of interface between the 
optical extent of the galaxy and the surrounding medium, 
is clearly important. Edge-on spiral galaxies, such as the 
galaxy we study here, NGC 891, are oriented in such a way 
that these are particularly useful in the study of this in- 
terface region dStrickland et alll2004D . Because of the en- 
ergies associated with supernovae and the wind of massive 
stars, observations at X-ray energies are of particular inter- 
est in studying the hot haloes or outflows from galaxies. Of 
the two currently operating major X-ray missions, Chan- 
dra has very high spatial resolution and moderate collecting 
area, whereas XMM-Newton has reasonable spatial resolu- 
tion, but a larger collecting area, and thus should be able to 
see lower X-ray surface brightness features. 

In this paper we will present results from XMM- 
Newton and Chandra observations of NGC 891, an edge- 
on spiral galaxy. NGC 891 is very similar in many respects 
to our own galaxy. The high inclination makes NGC 891 
ideal for studying the diffuse X-ray emission from the 
outflow extending from the plane of the galaxy (partic- 
ularly in the NW direction), which has previously been 
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studied with the ROSAT feregman fc Pildidll992Th ASCA 
jBreeman fe Pildisll99llRead Ponman fc Stricklandll997l) 
and Chandra ^S^iM^id^^Ij^On^r ol^servsLtones. 

High-resolution X-ray studies reveal diffuse emission up 
to about 10 37 erg s _1 , but also reveal a rich point source 
associated with stellar sources, such as X-ray binaries and 
supernovae. The X-ray point source population is closely re- 
lated with the hist ory of star formation in the galaxy (e.g., 
IKilgard et all2002lh The slope of the X -ray luminosity func- 
tion (XLF) of point sources is believed to be a good indica- 
tor of this, with flatter XLF slopes indicating more recent 
star-formation (due to a large fraction of high luminosity 
massive X-ray binaries). The point source populations of 
sever a l nearby galaxies have been analysed llHartwell et alJ 
[20041: IKilgard et alJ 120021 : IColbert et ail l2004h . Of them. 
lHartwell et alJ <l2004h selected a small group of well-studied 
spirals, and found that the XLF slope is generally steeper 
for normal spiral galaxies than for the starburst galaxies, 
and correlated with other observable parameters. In this 
work, we compare these X-ray, optical and IR properties 
of NGC 891 with a similar sample of normal and starburst 
galaxies. 

As part of the point source population, the presence of 
an active galactic nuclei (AGN) is also the subject of much 
debate in NGC 891 and other galaxies, and how dominant 
or irrelevant the A GN is to the energetics of the central re- 
gions of the galaxy. IStrickland et alJ i2004Tl show that there 
is the possibility of a weak hard X-ray source (2.0-8.0 keV) 
using the Chandra data. The discovery of a nuclear AGN 
in NGC 891 would help us to categorise it as a starburst 
galaxy. 

A rare form X-ray point sources in galaxies are X- 
ray luminous young supernovae. There have been few X- 
ray detections of young supernovae, and the ones that have 
been detected are associated with relatively nearby Type 
II SN events llHouck et alJ 1^998) and references within). 
One of these rare objects, detected at X-ray wavelengths, 
is SN1986J, located in NGC 891 and detected as a bright 
source in our X-ray observations. It is also very radio bright. 
iHouck et ai] Jl998fl studied the evolution of SN1986J using 
both ASCA and ROSAT data. SN1986J is still visible at X- 
ray energies and we present a new analysis of the evolution 
of this source using both Chandra and XMM-Newton data. 

The paper is organised as follows: In §2 we discuss the 
XMM-Newton and Chandra data used here, and present 
the reduction procedure and results of the analysis of point 
sources (Chandra) and the diffuse emission (XMM) in §3. 
In §4, a comparison of X-ray and near and far-IR properties 
is made with other nearby spiral galaxies. Observations of 
the supernova SN1986J are presented in §5, and compared 
with observations at previous epochs. General conclusions 
are presented in the final section. 



2 OBSERVATIONS 

NGC 891 is similar to the Milky Way in optical luminos- 
ity (m ^ = 9.4 from RC3, implying M B = -20.4), Hubble 
type (Sb) and rotat ional velocity (225 km s~ 1 . lvan der Kruhl 
119841: lRupenlll99ll) . It is a nearby edge-on spiral, and thus 
has been the subject of many detailed studies of interstel- 
lar dust and gas, from observations of the radio contin- 



Table 1. The nearby edge-on spiral galaxy NGC 891. 



Parameter Value 

RA (h m s) 02:22:33.41 

DEC (° ' ") +40:20:56.9 

Diameter (major) 14.1' 

Diameter (minor) 3.1' 

Major-axis position angle 22° 

Redshift 528 km/s 

Adopted Distance 9.08 ± 0.45 Mpc 



Notes: RA and Dec are given in J2000. The diameters 
correspond the D25 ellipse (i.e. the length of the axis at the 
isophotal level of 25 mag/arcsec 2 in the B Band) Source: NED, 
LEDAS. Source of distance explained in §2. 

uum (e.g.. iDahlem. Dettmar. fc Hummel 1 1994h . H I (e.g., 
ISwaters 1 Sancisi 1 fc van der Hulst 1997), carbon monoxide 
fe.g.. | So fue fc Nakail Il 993ll an d molecular hydrogen (e.g., 
IValentiin fc van der Werdll999lh 

However, there is considerably more star formation 
in NGC 891 than in the Milky Way, presumably due 
to the pr esence of abo ut 2.5 times as much molec- 
ular gas dScoville et all Il993l) . As a result, NGC 891 
is found to be twice as luminous at IR wavelengths 
iWainscoat. de Jong fc Wesseliusl H987). and there is evi- 
dence of enhanced d ust extinction for opt ical light, even 
outside the plane (e.g.. lHowk fc Savageil997l) . The extended 
gaseous halo of the galaxy in evident from Ha observations 
using HST/WFPC2, wit h filaments reachin g up to 2.2 kpc 
above the galactic plane iRossa et alJ l2004j. 

Other wavelengths have yielded other interesting de- 
tails. Pol arised radio emission ha s been detected from 
NGC 891 (ISukumar and Allenll99ih . illustrating the nature 
of the interstellar magnetic field. A luminous radio halo 
has been well established with a scale height of 2.7 kpc 
jAllen. Baldwin fc Sancisilll978l) . 

Throughout this analysis, we assume a distance of 
9.08 ± 0.45 Mpc for NGC 891, which is a weighted 
mean of the values obtained bvlCiardullo. Jacobv fc Harris I 
Jl99ll) and iTonrv et all toOlR Other selected parame- 
ters are shown in Table Where relevant, we use Hq = 
72 km s -1 Mpc -1 . 

2.1 XMM-Newton observations 

NGC 891 was observed with the XMM-Newton telescope on 
August 22 2002, as part of a GTO observation. The galaxy 
was at the centre of the field of view, with the entire emission 
also contained within the field of view. 

The XMM-Newton data was observed with three dif- 
ferent cameras, MOS1, MOS2 and PN. Each of the data 
sets were analysed separately. After running the standard 
SAS tasks: cifbuild and odfingest; the chains to process 
the odjs were run (emchain, epchain) simultaneously run- 
ning the badpixfind command. These have proved to be 
considerably better at detecting bad pixels than the tasks 
EMPROC and EPPROC. Having extracted the lightcurve in the 
10-15 keV energy range, the observation was found t contain 
a considerable amount of flaring which we used to identify 
the good time intervals. The periods of high background 
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Figure 1. a): The XMM-Newton soft X-ray (0.3-2.0 keV) contours overlaid onto the optical Digitised Sky Survey (DSS) image. The 
lowest contour level is set at 3a above the background, b): point sources. 



were identified by eye at 55/40.0 counts for the MOS/PN 
cameras respectively. The resulting good times are 14.6 ks 
(originally 18.0 ks) for MOS1/MOS2; and 9.6 ks (originally 
15.0 ks) for PN. 

Events were filtered based on the the pattern parameter, 
which indicates the geometry of the detection of each event, 
i.e. the number of adjacent pixels that detect each photon. 
The events that were retained in the filtering were the ones 
that are well-calibrated (singles, doubles and quadruples for 
MOS & singles and doubles for PN). Due to the extent of 
the diffuse emission present in this galaxy, we used blank 
sky background files for subtr action of our data. We fol- 
lowed the technique set out bv lRead fc Ponmanl (l2003fi for 
producing images and spectra that are correctly background 
subtracted. 

Given the high density of point sources, and the high 
inclination, it seemed that manually setting a region with 
exclusion regions based around the point sources (extracted 
from Chandra data- see below) was more appropriate. The 
radius of the exclusion region was set by the SAS task 
CALVIEW, using the encircled energy (PSF) function, at 15" 
(corresponding to a PSF fraction of 75%). When the sources 
were removed, we used the CIAO task dmfilth to interpo- 
late over the holes. The width of the background annuli were 
initially taken to be 1.5 times the radius of the source region, 
but many were adjusted to ensure that the background an- 
nuli did not overlap with any source regions. In some cases it 
was impossible to extract adequate background annuli due 
to the compact nature of the sources on each other. In that 
case, the sources were merged and treated as a single source 
removal region. When analysing the spectral data, the data 
from the three cameras were fitted simultaneously together 
to ensure the best fit. However, for the radial profile analysis, 
the data was mosaicked together using the task EMOSAIC. 



2.2 Chandra observations 

To better characterise the point sources, we used the 
archived Chandra/ 'ACIS-S observation (51 ks, PI: Breg- 
man), which was reduced in a similar way, using the CIAO 
software with online threads. Once the Chandra data prepa- 
ration was completed, the CIAO tool wavdetect was run 
on the Chandra data. We narrowed down the returned list 
of sources by excluding all the sources outside the D25 el- 
lipse (see Table 1), and considering only the sources inside. 
We used the wavelet scales 1, 2 and 4 to extract the point 
sources, as these provided the best identification of all the 
point sources present. These sources are shown on a DSS 
overlay in Fig. 0. We used the results to generate an X-ray 
flux of each source to generate a luminosity function (Fig.|3J. 



3 ANALYSIS OF THE X-RAY OBSERVATIONS 

Here we analyse the X-ray observations described in the pre- 
vious section to characterise the population of point sources 
in NGC891, using the archived Chandra data set, utilising 
the superior resolution of Ciandra/ACIS setup. We use this 
observation to subtract the point sources from the XMM- 
Newton observation, to investigate the nature of the diffuse 
emission. 



3.1 The Point Source Luminosity Function from 
Chandra 

We detected a total of 26 point sources within the D25 el- 
lipse, down to a flux limit of 10 -15 erg cm -2 s _1 , which 
at the adopted distance of NGC 891 amounts to a lumi- 
nosity of F m i n = 10 37 erg s" 1 . The regions were all output 
as ellipses containing 99.7% (3er) of the source counts. As a 
considerable number of the sources contained few counts, an 
accurate spectral fit was not possible, so we modelled each 
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Figure 2. The X-ray luminosity function (XLF) for point sources, in NGC 891, in the 0.3-8keV (left) and 2—8 keV (right) energy ranges, 
measured from the Chandra observation. The slope, defined in 0, is calculated to be a = 0.77^q '\q and a = 0.77±g;^ respectively. 
Since the slope does not substantially change for the higher energy range, we conclude that differential extinction within the inclined 
spiral does not significantly affect the shape of the XLF. The inset in each plot shows the shape of the log likelihood function, derived 
from the method of evaluating the slope and errors, discussed in §3.1. 



of the point sources individually in XSPEC using an ab- 
sorbed power-law model with a slope of T = 1.8. The point 
sources with more than 100 counts were individually fitted. 
Flux values were established for each of the point sources in 
an energy range of 0.3-8.0 keV, and these were combined to 
make a luminosity function (Fig. |2Jl. We discuss the point 
source associated with SN1986J in greater detail in the fol- 
lowing section. 

As we will see below, the slope of the luminosity func- 
tion can be used as an indicator of recent star formation 
activity in a galaxy. This so-called log JV-log S relation is 
usually expressed in the cumulative form 



log N(> S) = -a log 5" + k 



(1) 



where a straight line is fitted to the cumulative his- 
togram. Here we adopt a more robust method of mea- 
suring the slope of the luminosity function with the po- 
tential of better measuring the uncertainty of the slope. 
This is a gen eralised version of the much-u s ed m ethod 
adapted from ICrawford. .Tauncey. fc Murdoch! l)l970h and 
iMurdoch. Crawford, fc Jauncevl lll973t) . the crucial differ- 
ence being that in the former, errors in flux measures are 
not taken into account and in the latter, all errors are con- 
sidered equal. 

Here, each of our n point sources has a measured lumi- 
nosity Fi, given the adopted distance to the galaxy, and an 
independently estimated error a%. We represent the proba- 
bility distribution of the luminosity S of point sources in the 
galaxy in its differential form 

P(S) dS = AS' 13 dS. (2) 

On comparison with Jy, a = /3 — 1. Our exercise thus con- 
sists of maximising the log likelihood function 



(3) 



where the distribution of our measured values of flux and 
standard deviation (Fi,t7i) is given by 



P(Fi,ai) = 



f™P(Fj,(Ti\S)P(S) dS 
f™P(F i ,<T i \S)P(S)dSdF' 



(4) 



Assuming the errors of measuring flux and luminosity are 
distributed as a Gaussian, the integrand above is given by 



P{Fi,a t \S)P(S)dS = 



S~ p e 



!3-{F i -S) 2 /2a 



dS. 



(•>) 



We numerically find the value of (3 (thus a) for which 
C in © is maximum. For large N, the probability distri- 
bution for p is asymptotically Gaussian, and the lcr error 
in [3 corresponds to AC = 0.5. For our limiting 0.3-8 keV 
luminosity, F m i„ = 10 37 erg s _1 , we find the slope of the 
cumulative LF given by CJ to be q = 0.77t ;", as shown 
in Fig. H 

Since NGC 891 is almost exactly an edge-on galaxy, 
the effect of differential local extinction might be impor- 
tant within the point source population, which might lead 
us to progressively miss the fainter objects on the far side 
of the galaxy. Since the effect of such extinction is energy- 
dependent (being larger at lower energies, particularly be- 
low 2 keV), we tested whether it is substantial by re- 
calculating the luminosity function for the energy range 2- 
8 keV (Fig. |2t>)- The slope was found to be the same at 
q = 0.77t u 'i2, leading us to conclude that the effect of dif- 
ferential extinction is not very serious for this target. 



3.2 The diffuse X-ray emission from 
XMM-Newton 

We extracted the spectrum of the diffuse emission in the 
energy range 0.3-6 keV, excluding a circle of 15 arcsec 
around each detected point source, from the 3 XMM-Newton 
camera data and analysed them with the heasoft package 
XSPEC. We performed simultaneous fitting of all three data 
sets in order to fit the best model. We also tried to add the 
Chandra diffuse emission to the fits, but due to the lim- 
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Figure 3. (a): The X-ray spectra of NGC891, with an absorbed single temperature model + power law. Spectra for all three XMM- 
Newton cameras have been simultaneously fitted in the 0.3-6.0 keV energy range. The top spectrum is from the PN camera, (b): As in 
(a), but the fitted model is an absorbed dual temperature + power law. 



Table 2. Summary of recent fits to the diffuse X-ray emission. 
The XMM-Newton data are the results generated from this work. 



NGC891 minor axis surface brightness profile 



Obs. 


Model 


feTi 


kT 2 


Y 2 






(keV) 


(keV) 


XMM 


pow+mek 


0.26 ± 0.01 




1.06 


XMM 


pow+mek+mek 


0.08 ±0.01 


0.30 ±0.03 


1.05 


ROSAT"- 


mek+mek 


0.31 


10 




ROSAT* 


mek 


0.11 ±0.03 




1.10 



Notes: References: a) IBregman fc Pildisl Il994ft . b) 
Read. Ponman & Stricklan dl lll997D 



ited collecting area of the Chandra telescope, we found the 
diffuse emission spectra generated did not improve the fit. 

We initially fitted an absorbed single temperature 
plasma model with a power-law component to the data. We 
fixed the Galactic value of column density rig of 6.78xl0 20 
cm -2 (obtained from Golden). Since the target is an edge- 
on galaxy, the local column density is expected to be 
substantial- this was looked for, and fixed at a fairly low 
value of tih (local) at 2.0xl0 19 cm -2 . The power-law model 
was added to our data to account for point source contami- 
nation. 

Figure^ shows the best fit to all three data sets, hav- 
ing a x 2 =68.35 with 65 d.o.f., leading to Xred=l-06. We 
extract a value of kT = 0.26 ± .01 keV. In comparison, 
iRead. Ponman fc Strickland! <ll997l) fit a single temperature 
model to the dat a with a cooler tempera ture component 
(0.11 ± 0.03 keVl. lBregman fc Pildisl lll994t) found that the 
data was best fitted by a two-temperature model with kT 
at 0.31 and 10 keV. 

As the result from IBregman fc Pildisl (119941) implied a 
two temperature model, a second temperature component 
was added to the data (Fig. |3>). The two temperatures 
that were inferred, as a result, were kT = 0.08 ± 0.01 keV 
and 0.30 ± 0.03 keV. The two temperatures agree within 
errors to the coo l er tem perature component derived by 
IBregman fc Pildisl <|l99 4|) and the single temperature com- 
ponent derived bv lRead. Ponman fc Strickland! dl997t) . 




-5000 5000 

Distance from centre (parsec) 



Figure 4. Minor axis radial surface brightness profile of the dif- 
fuse emission. The plot is in the 0.3-8.0 keV energy band. 



However, it should be stated that the spectrum of the 
data was best constrained by a single temperature fit, and 
the improvement in the Xred by adding a second tempera- 
ture component is not statistically significant (f-test statistic 
value = 1.05). 

A summary of our fit param e ters com- 
pared with those of IBregman fc Pildisl (|T994) and 
IRead. Ponman fc Strickland! i^9\^i are presented in 
Table H 

Our results do seem to complement the ROSAT data to 
some extent. Our single and dual temperature model values 
are consistent with temp eratures that both have derived. 
IBregman fc Pildisl Jl994li do find a high temperature com- 
ponent at 10 keV, which is not determined from our fits. The 
reason for this may be the poor resolution of ROSAT, which 
means that it is difficult to remove point sources effectively. 
This effect is compounded by the high inclination of the sys- 
tem. Also, SN1986J was considerably brighter, at the time of 
the ROSAT observation more than ten years ago, than it is 
now, and may well have contaminated the extracted diffuse 
emission. Finally, the hard X-ray flux of a nearby unresolved 
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bright point source near S N1986J might al so acc ount for this 
discrepancy (iRead. Ponman fc Stricklandlll997|) . 

We also extracted images of the diffuse emission from 
the XMM-Newton data. These i mages complement the X - 
ray observati ons from ROSAT feregman fc Pildisl Il994l) 
and Chandra iRead. Ponman fcStoicklandl^b^ . There is 
a considerable amount of emission from the centre of the 
galaxy in the NW direction. Figure^ shows this result for 
the soft energy band (0.3-2.0 keV). The X-ray contours are 
overlaid onto the optical Digitised Sky Survey (DSS) im- 
age. The image was smoothed with an adaptive gaussian 
program ASMOOTH such that the ler width of the gaussian 
contained 40/3 counts around each pixel. We also took a 
minor axis diffuse emission profile to show the extent over 
which the extended emission protrudes (Fig. ■ The emis- 
sion protrudes out to just over 6 kpc from the centre of the 
galaxy in the NW direction, but drops off rapidly in the SE 
direction outside the plane of the galaxy. 



3.3 Evidence of central AGN 



Using the Chandra observation, IStrickland et al.l (120041) 
found a weak hard (2.0-8.0 keV) X-ray source near the nu- 
cleus of the galaxy, de fined as the p osition of the radio con- 
tinuum point source iRupenl ll991). This was suggested to 
be a weak AGN possibly associated with the radio source. 
Unfortunately, we were not able to detect this source with 
our XMM-Newton observation, most likely due to the poorer 
resolution compared to Chandra. 



4 COMPARISON WITH OTHER LOCAL 
GALAXIES 

In order to characterise the global properties of NGC 891 
and find whether its properties are consistent with that of 
a starburst galaxy or not, we compile various X-ray, near- 
IR and far-IR properties of several nearby spiral galaxies, of 
which some are normal galaxies and some are actively star- 
forming and starburs t galaxies. Star t ing fr om a sample of 
gala xies selected from [ lartwell et al\ (| 2004|). a few gala xies 
from IStrickland et al.1 <l2004h and lHartwell et al.1 (l2004h . A 
selected list of parameters are collected in Table |S| 



Far-Infrared luminosities We calculated the FIR lu- 
miosities (erg s _1 ) differently depending on whether the 
source was fro m ISO (preferably) or IRAS. We used the 
equations from Ben do et alJ |2002) (fluxes in Jy): 



L FIR = 1.89 x 10" 
and for IRAS 



(1.36/eo + 0.958/iQo + 0.439/iso), (6) 



L FIR = 1.26 x 10- 14 (2.58/ 60 + /ioo), 



(7) 



Near-Infrared luminosities To calculate the K- 
band luminosity Lk (erg s _1 ), we derived the fol- 
lowing expression using the ze ro point values from 
ICohen. Wheaton fc Meeeathl <2003ft : 



log L K = 43.13 - 0AK tot + 21ogL>, 



(8) 



with D in Mpc, and K to t the t otal K-band m agnitude taken 
from the 2MASS galaxy atlas Jjarrett et alJl2003l) . 



X-ray luminosities Lx was established from the liter- 
ature, but different au thors used varying energy ranges. 
IStrickland et alJ i2004f) evaluate their fluxes for an energy 
range 0.3-2.0 keV (soft X-ray) We have adopted this en- 
ergy range to keep our values for Lx ■ In the case where the 
data was not in the necessary energy range, we used the 
Portable Interactive Multi-Mission Simulator (PIMMS) to 
generate our data. Using the model fit parameters for the 
diffuse emission in the literature, we estimated a revised flux 
in the appropriate energy range. 

Deprojected area of galaxyWe differ from the calcula- 
tion made_bvjHar^eU_eVaL| i2004l) for the area of the D25 
ellipse. lHartwell et al.l l|2004|) calculated the galaxy area as 
nab cos#. For systems that are face on (such as NGC 4214) 
this calculation is adequate. However, for very high incli- 
nation systems, the measure of the semi-minor axis simply 
represents the width of the disk, and cannot be used in cal- 
culating the area. Therefore, for the sake of consistency we 
choose to evaluate our area as ira 2 for all galaxies in our 
sample. 

4.1 Starburst or not? 

There has been much discussion as to whether NGC 891 
is a normal spiral ga laxy, or whether it is a starburst 
JStrickland et al.ll2004b . We try to categorise NGC 891 by 
means of comparison with a sample of 16 nearby galaxies. 
Figure contains four plots, which are described below. 

The slope of the luminosity function can be used as 
an indicator of the general properties of the host galaxy. 
If the star formation rate in the galaxy is constant, then 
the luminosity function of the sources should be fit by a 
single unbroken power-law model. If the galaxy is a starburst 
one, then new high mass X-ray binaries (HMXBs) would be 
formed, breaking the luminosity function slope. The break 
in the slope would decrease with time, and would be an 
indication of the time of previous bursts in the galaxy. 

As we saw above, in NGC 891, we chose to fit the XLF 
slope can be fitted by a single power-law. In the case of 
NGC 1482, the galaxy is further away (22 Mpc) and only two 
sources could be detected so no XLF slope was determined. 
Two other galaxies NGC 4244 and NGC 6503 also have low 
statistics (with 3 an d 4 sources respective ly), but XLF slopes 
were calculated bv lColbert et all i2004l) . We include these 
in our plots, but they do not appear to be representative of 
the normal spirals category. 

Figure plots the XLF a gainst the Seo/Sioo flux ra- 
tio. [Lehner^&^eckman] lll996l) analyse the properties of su- 
perwinds and quantitatively discuss the likelihood of gener- 
ating a superwind with respect to the infra-red properties 
of the galaxy. The galaxies with high luminosity (Lir ^ 
10 44 erg s _1 ), large infra-red excesses (Lir/ 'Lopt ^ 2) and 
warm far infra-red colours (Seo/Sioo ^ 0.5) are most likely 
to have galactic superwinds. For the sample selected, the last 
criterion is satisfied for all the starburst galaxies. Of the 
starburs t galaxies in our samp le, superwi nds are detected 
in M 82 (IStrickland et alJll997l). NGC 253 (IStrickland et"al] 
l200d) and NGC 4945 jMoorwood et alJ Il99rj). NGC 5253 
JStrickland fc Stevensll999D and NGC 4449 JSummers et alJ 
l2003Tr exhibit emission from superbubbles. Although super- 
winds and superbubbles are both the result of massive star 
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Figure 5. Top left a): the XLF slope plotted against the IRAS flux ratio S60/S100 for a sample of nearby galaxies. The circles represent 
starburst galaxies; crosses represent normal spirals; square for interacting galaxies; and the filled triangle represents NGC891. The same 
symbols are used throughout this paper. The data labels have been removed for all other plots for neatness. The only labels left in are 
those referred to in the text. The error on NGC 891 is the standard deviation, discussed in §3.1. Top right b): the XLF plotted against 
far infra-red luminosity Lfjr. Bottom left c): the XLF plotted against the ratio of the IRAS 60/i flux and the spiral galaxy area. Bottom 
right d): the IRAS flux ratio plotted against Lpjn normalised with the K-band luminosity. 



formation processes within the densest regions of the hot 
galaxies, the winds are able to channel the metals produced 
straight into the IGM, whereas the superbubbles do not 
reach the outskirts of the host galaxies, hence retaining their 
newly process ed metals, consequently raising the abundanc e 
of the ISM ijTenorio- Tade. Silich fc Munoz-Tu ianl 1200^ . 
The only starburst galaxies not to have superwinds or super- 
bubbles are NGC 4214, which is calculated likely to exhibit 
blowout l|Hartwell et, al.ir2004l) and the Antennae, which is 
an interacting system, so no superwind has formed due to 
disruption from the merging. 

Therefore, there are two distinct categories of result: 
The low XLF slope galaxies have a higher ratio of warm IR 
colours than the high XLF slope galaxies, which are normal 
spiral galaxies. NGC 891 appears to have the XLF slope of 
a starburst, but the S60/S100 corresponding more to normal 
spirals than starbursts. 

Figure QJ) plots the far infra-red luminosity (Lfir) 



against the XLF slope of the data. Lfir is an indicator 
for how much dust there is in a galaxy, and hence the star 
formation rate teennicuttJl2003l) . The FIR also contains the 
most important cooling lines in the neutral ISM. Therefore, 
if the XLF slope is dependent on star formation rate, then, 
one would expect a correlation between these two param- 
eters. However, this is not observed in our sample. Both 
the starbursts and normal galaxies span a wide overlapping 
range of FIR luminosities. The division between spirals and 
starbursts is quite clear in this plot, breaking either side of 
1. Therefore, an XLF-Lfir dependence cannot be justified 
on the basis of these results. 

Figure^: plots the XLF slope against the S60 nu x scaled 
against the area of the D25 ellipse. The x-axis therefore repre- 
sents the star formation rate per unit area. There is a definite 
correlation between the data points on this plot. The higher 
star formation rate corresponds to starburst galaxies, which 
we would expect to see. The major exceptions to this trend 
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are NGC 4244 and M83. The XLF slope of NGC 4244 was de - 
termined with only three data points JColbert et al]l2004) . 
and has an estimated error on the data. M83 should be clas- 
sified as a starburst, as observations indicate t hat there are 
mass ive clusters of stars in th e nuclear region ([Har ris et alJ 
2001). However, as discussed in iKileard et aH (2002) the star 
formation rate is low compared to other starbursts, and the 
starburst region is confined to a small area of the galaxy. 
Therefore, most of the point sources analysed were taken 
from a region outside of the starburst region, so the XLF 
slope is more representative of the disk population of nor- 
mal spirals as opposed to the standard starburst population. 

In Fig. 01, the warm IR colour ratio S60/S100 is plotted 
against the ratio of Lfir/ Lk ■ This ratio tells us the extent 
of the star formation activity normalised against the stellar 
types present in the galaxies. The different galaxy types are 
clearly grouped separately, and all the galaxies seem to fol- 
low a linear trend. The normal spiral galaxies have lower star 
formation rates, which would imply that they are less likely 
to generate galactic superwinds. The galaxies with higher 
star formation rates are all starburst galaxies with super- 
winds (or superbubbles) present. Again, NGC 891 seems to 
sit between the two categories, making it difficult to cate- 
gorise as either a spiral or as a starburst. 

4.2 The X-ray Schmidt law 

In a classic paper, ISchmidil il959l) showed that the rate of 
Population I star formation is proportional to the density 
of available gas in a galaxy. This has been re-discovered in 
various guises over the last few decades. In terms of the X- 
ray luminosity of the galaxy, which is related to a density of 
hot gas, one expects it to correlate with the star formation 
rate , which we can characterise by th e total far-IR flux. 

iRanalli. Comastri fc Settil (|2003j) analyse the effective- 
ness of the 2-10 keV X-ray luminosity as a star formation 
rate indicator, showing the linear relation between the X-ray 
luminosity and the far-IR luminosities, although there are 
discrepancies between many other published results. Previ- 
ous studies using ASCA, ROSAT/PSPC or EINSTEIN val- 
ues would overestimate Lx, being unable to separate out 
the emission from X-ray binaries and AGN. 

iGriffiths fc Padovanil il99Cf> established L x oc L° F 6 IR , 
although 560 made up the x-axis, and Lx ranged 
from 0.5-3.0 keV. Als o with the EINSTEIN satellite, 
IDavid. Jones fc Formanl jl992l) obtained an almost linear 
fit to their data Lx oc Lp'j'^ ' 06 agreeing quite well with 
our results within errors. This result was determined in the 
0.5-4.0 keV energy band and t he sample was a selection 
of star burst and normal galaxies. IRanalli. Comastri fc Settil 
( 2003) analyses data obtained from the ASCA satellite 
and for an soft X-ray energy range of 0.5-2.0 keV obtains 

n ft*7— i— n 

Lx oc L FIR ' , an identical result to ours for a very similar 
energy range (we use 0.3-2.0 keV). 

We plot Lx (the diffuse X-ray luminosity) against Lfir 
for our sample galaxies in Fig. |S] We establish that the re- 
lation between the data is Lx oc Lp^ 7 ^ ' 07 , without the 
two data points on the plot which do not appear to lie near 
the line. If these are included, the fit becomes considerably 
worse: Lx oc Lp 71 ^ 10 ' 12 . These are the two starburst galaxies 
NGC 253 and NGC 4945, both of which contain a superwind 
emitting from the centre. This would indicate the presence 
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Figure 6. The X-ray Luminosity against IR luminosity for a sam- 
ple of galaxies. The circles represent starburst galaxies; crosses 
represent normal spirals; the square for interacting spirals; and 
the filled triangle represents NGC 891. 

of an excess of dust in these galaxies containing considerable 
star formation. 

NGC 891 lies comfortably near the best-fitting line, 
away from the normal galaxies since its X-ray and far-IR 
fluxes are higher, again showing that it does not have as 
extreme a rate of star formation as a starburst galaxy. 



5 THE RESIDENT YOUNG X-RAY 
SUPERNOVA SN1986J 

Even though thousands of supernovae have been found 
in the optical, there are instanc es of only 15 super- 
novae being detected in the X-ray ilmmler fc Lewinll2003l : 
IBreeman etai] i2003). Of these, a handful have been mon- 
itored over many years, their X-ray emission dominating 
their radiative outp ut after about a year after explosion 
iPoolev et alJ 12002'). Thermal emission from the reverse 
shock region is expected to be seen as softer X-ray emis- 
sion within the expanding SN shell a s it interacts with the 
dense stellar wind of the progenitor dChevalier fc Franssonl 
11994) . 

One of these Type II supernovae observed in the X- 
rays is SN1986J in NGC 891, whi ch was discovered using 
radio observations from the VLA (IRuoen et alJll987f) . The 
first X-ray observations o f SN1986J were taken with ROSAT 
feregman fc Pildisll992^ . Since then, two more ROSAT and 
two AS CA observations h ave been undertaken in the mid 
1990's jHouck et al.lll998lh The X-rays are thought to be 
caused by a shock pro pagating through the former stellar en- 
velope dBregman fc Pildislll992h : or due to a the supernova 
envelope colliding with a shocke d clumpy wind, whi c h pro - 
duces the X-rays <ChugaJll993h . iBregman fc Pildisl (Il992ft 
found that the X-rays were emitted primarily in the soft 
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Figure 7. The Chandra, X-ray spectrum for SN1986J. 



X-ray band (0.1-2.5 keV). An an alysis of the ROS A T and 
ASCA observations was done bv IHouck et alJ (^998), who 
determined that the X-ray light curve of the data was drop- 
ping by t~ 2 , which was considered to be a rather fast decline. 

A summary of all observations are shown in Table [3] 
Note that the X-ray fluxes derived for the ASCA data are 
an average over all their different models, and so the error 
shown here is simply the standard deviation of these data 
sets. 

5.1 The spectrum of SN1986J 

The spectrum for the Chandra observation was extracted us- 
ing a local background and we fitted different temperature 
models to the data. Firstly, an absorbed single temperature 
thermal plasma model (mek al) was fitted with a loc al and 
Galactic absorbing column fereeman fc Pildisl 17992) . The 
Galactic value for lognn was fixed at 20.8 (tih in cm -2 ) and 
the local absorbing column was allowed to vary. The abun- 
dance was initially fixed at 0.3 then 1.0 solar, and finally was 
allowed to vary. The local logn.tr was evaluated at 21.6, in 
good agreement with the value of iBreeman fc Pildisl <ll992l) 
(21.7), but the temperature range that we obtained from 
the model was 3.5-4.1 keV. This temperature component 
is much larger than the 2.0 keV values obtained from the 
ROSAT observations, but smaller than the ASCA tempera- 
ture ranges. A vmekal model was also tried, where the indi- 
vidual element abundances are allowed to vary, as the mekal 
fit results were leaving several emission lines unidentified. 
The result from the vmekal fit is plotted in Fig. Q The tem- 
perature was determined to be 4.6 keV, local log tlh =21.5, 
and X-ray flux 0A5to'.l x 10~ 14 erg s _1 cm~ 2 . Many elements 
were identified in the spectrum of the supernova, most no- 
tably, the Fe emission line at 6.4 keV. The abundances are 
tabulated in Table 21 

On inspection of Fig. Q some of the emission lines are 
not fitted that well by a vmekal model. Several methods 
were tried to improve the fit, including a two temperature 
model (mekal+vmekal) and a different single temperature 
thermal plasma model (apec). In each of these models, the 
hard temperature and local nn were similar to the results 
obtained from the single temperature fitting. For the two 
temperature fitting, the lower temperature component was 
0.11 keV which is likely the component from the diffuse emis- 



Table 3. Summary of X-ray observations of SN1986J. 



Telescope Date MJD Flux 

(1) (2) (3) (4) 



ROSAT 


Aug 1991 


48486.1 


7.90 ±0.41 


ROSAT 


Jul 1993 


49206.5 


4.91 ± 0.24 


ASCA 


Jan 1994 


49374.5 


6.07 ±0.13 


ROSAT 


Jan 1995 


49744.0 


3.95 ± 0.26 


ASCA 


Jan 1996 


50113.5 


5.28 ±0.10 


Chandra 


Nov 2000 


51852.0 


96+ 05 


XMM-Newton 


Aug 2002 


52508.0 


0.85 ± 0.15 



NOTES: Column (1): The first 2 ROSAT observations were 
taken with the PSPC, the last with the HRI. Column (4): The 
X-ray fluxes are measured in 10 — 13 erg s —1 cm" 2 . The ASCA 
and ROSAT data are compiled from lHouck et all ll998l) : the 
Chandra and XMM-Newton fluxes are from this paper. 



Table 4. Metal abundances for single VMEKAL model. 



Metal 


Abundance 


Ne 


5.12 


Mg 


2.71 


Al 


5.42 


Si 


1.87 


S 


3.8 


Ar 


1.09 


Ca 


3.48 


Fc 


0.92 



NOTES: Abundances are relative to solar. 

sion of the galaxy, and not the supernova. The results from 
all the Chandra fits are tabulated in TableQ)] None of the fits 
identified any of the other emission lines. As all the flux val- 
ues generated from the fits were similar, the value adopted 
for our supernova flux is the weighted mean of the different 
models: 9.6lg \ x 10 14 erg s _1 cm -2 . The energy range cho- 
sen was 0.5-2.5 keV to be consistent with the energy ranges 
used in the ASCA and ROSAT data. 

For the XMM-Newton data, the flux was much harder 
to determine. As the flux is dropping every year, the su- 
pernova becomes harder to model the later the observations 
are. We fitted all three cameras simultaneously to try and 
obtain the best fit to the data. Our data was best fitted with 
a single temperature model (mekal) with a temperature of 
3.6 keV and absorbing column density lognn = 21.53. This 
fit resulted in a flux of 8.31 ± 1.06 x 10 14 erg s" 1 cm~ 2 . 
Fitting a dual temperature model did not improve the fit. 
The temperatures obtained from the fit were 0.21 keV and 
3.82 keV for the soft and hard components respectively. This 
fit did produce a surprisingly large local uh of 22.6. The two- 
temperature flux was 8.40±1.16 x 10~ 14 erg s^'cm" 2 . The 
value adopted from the flux is the mean of the four differ- 
ent temperature fits for the XMM-Newton data: 8.5±0.15 x 
10" 14 erg s _1 cm" 2 . 

5.2 The evolution of SN1986J with time 

SN1986J has been observed with ROSAT, ASCA, Chandra 
and XMM over the last twelve years. In Fig. |H|we plot the 



10 R. F. Temple, S. Raychaudhury and I. R. Stevens 



Table 5. XMM and Chandra single-temperature fits to the su- 
pernova lightcurve. 



Telescope 


Model 


kT 


nH 


Flux 


x 2 rf 


Abund 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Chandra 


mek 


3.5 


0.41 


9.7 


1.34 


0.3F 


Chandra 


mck 


4.0 


0.39 


9.6 


1.16 


1.0F 


Chandra 


mek 


4.1 


0.39 


9.7 


1.15 


1.30 


Chandra 


vmek 


1.6 


0.30 


9.5 


1.09 




Chandra 


mek+vmek 


1.6 


0.34 


9.5 


1.13 




Chandra 


apec 


4.0 


0.39 


9.9 


1.39 


1.48 


Chandra 


vapec 


4.1 


0.37 


9.6 


1.18 




Chandra 


apec+vapec 


5.9 


0.45 


9.5 


1.16 




XMM 


mek 


3.6 


0.35 


8.3 


0.92 


0.77 


XMM 


apec 


3.4 


0.37 


8.6 


0.92 


0.77 


XMM 


vmck 


3.8 


1.15 


8.4 


0.90 




XMM 


mek+mek 


1.2 


0.38 


8.7 


0.84 





NOTES: Col (3): Temperature component of the supernova 
fitting in keV. For two-temperature models this is the harder 
component. Col (4): Absorbing column parameter local to 
NGC891 Xl0 22 cm- 2 Col (5): X-ray flux xl0~ 14 erg s^cm- 2 . 
Col (7): Abundance (relative to solar) from different fitting 
models. 'F' indicates a fixed fitting value. These are not shown 
for the vmekal/vapec model. 

0.5-2.5 keV flux against time over thi s period. We reca ll 
that from the ROSAT and ASCA data jHouck et al.lll998h . 
the decline of this flux was found to be oc t~ 2 , which was 
considered to be rather steep. 

Here we attempt to find a mean rate of decline over a 
much larger time scale. The data from ASCA and ROSAT 
appear to be inconsistent (iHouck et alJ Il99fetl and so we 
have established a best fit line inclusive and exclusive of 
the ASCA data. The slopes are calculated as —2.99 ± 0.45 
and —2.89 ± 0.19 respectively. In Fig. |5] we also show an 
alternative set of points for ASCA where we attempt to cor- 
rect for the poor resolution by subtracting from the the flux 
of the SN the combined flux of all the other point sources 
that are expected to be within its point-spread function. 

However, our results imply Lx oc t~ 3 . We can com- 
pare this with the expected behaviour from the work of 
IChevalier fc Franssonl 1^994), who show that Lx oc t^ 1 
when free-free emission dominates (T > 4 x 10 7 K), and 
Lx oc t~ ' 7 when line emission dominates (10 5 < T < 
4 x 10 7 K). Such a trend in the decline in luminosity is seen in 
various other Type II SN e, including SN 1993J i n M81 and 
SN1999em in NGC 1637 lllmmler fc LewirfeOOSt . However, 
our measured temperature of the plasma is around 4 keV, 
which corresponds to the former range, which indicates that 
the decline of Lx is unusually steep in this case. It is likely 
that the poor resolution of the ASCA and ROSAT satellites 
might not have been able to fully eliminate point source 
contamination, thus overestimating the fluxes of the data. 



6 CONCLUSIONS 

In this paper, we have presented a GTO XMM-Newton ob- 
servation of NGC 891, a nearby edge-on spiral galaxy. We 
find that the diffuse X-ray emission protrudes from the disk 
in the NW direction out to approximately 6 kpc. The ex- 
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Figure 8. The X-ray lightcurve of SN1986J. ROSAT data are 
shown by triangles; ASCA data by circles; Chandra data denoted 
by a cross; and XMM-Newton by a square. The ASCA data has 
been adjusted to remove point source contributions (see text). 
The lower values represent the point source corrected results. The 
x-axis represents the time in days since the supernova went off 
(t = 0, which is taken to be 1983.0). All data is in the soft energy 
range 0.5-2.5 keV. The solid line is the best fit line for all the data 
with adjusted ASCA data points. The dashed line is the best fit 
for the data excluding al l ASC A points. The dotted line is the 
best fit from lHouck et alJ l ll998ft for ROSAT only. Bottom panel: 
Temperature. 



traplanar hot gas was previously observed in ROSAT obser- 
vations, but the XMM observation seems to show a sharp 
cut-off to this gas, showing that this might be the extent 
to which the gas has reached. We also find that the best 
fitting thermal plasma model seems to require two temper- 
atures of 0.08 and 0.3 keV respectively, though the fit of 
a single-temperature plasma of 0.26 keV isn't much worse. 
This co ntradicts previous findings based on ROSAT ob ser- 
vations fereeman fc Pildislll99l iRead fc Ponmanll2003l) . 

We use an archived Chandra observation to study the 
point source population within the Z?2 5 ellipse of the galaxy. 
We use a robust maximum likelihood method to determine 
the slope of the cumulative luminosity function A^(> 5*) = 
S~ a , and find that the slope is rather shallow, a = 0.771q'jq. 
We have verified that this isn't predominantly due to ex- 
tinction, by plotting the XLF for the energy range > 2 keV, 
without much change on slope. 

Using a sample of other local galaxies, we have com- 
pared the X-ray and infrared properties of NGC 891 with 
those of nearby 'normal' and starburst spiral galaxies. We 
conclude that NGC 891 has more abundant star formation 
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than a normal spiral, but does not have as extreme proper- 
ties as starburst galaxies like NGC 253 and NGC 4945, and 
that it is most likely a starburst galaxy in a quiescent state. 

We examine an X-ray version of the "Schmidt Law", 
which correlates the rate of star formation in a galaxy to 
the mass or density of its available gas. We show that the 
diffuse X-ray luminosity, an indicator of gas mass, of nearby 
spirals scales with their far infra-red luminosity, as indicator 
of dust mass and star formation rate, as Lx oc L^;^ ±0 ' 07 , 
except for extreme starbursts, and NGC 891 does not fall in 
the latter category. 

We study the supernova SN1986J in NGC 891 in both 
XMM-JVewton and Chandra observations, nearly twenty 
years after its explosion. It was studied in the X-ray us- 
ing ROSAT and ASCA almost a decade ago. The flux and 
temperature calculated in those studies of low spatial reso- 
lution could have significant contamination from their point 
sources and the diffuse emission. A direct comparison shows 
that the temperature of the SN remnant at the time of ob- 
servation was 5 x 10 7 K, which is less than it was a decade 
ago, but it also reveals that the X-ray luminosity has been 
declining with time {Lx oc t~ 3 ) far more steeply than ex- 
pected. 
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Table 6. Luminosity function and flux data for a sample of nearby galaxies. 



Galaxy Distance LF slope 100 nm flux Source •Sgo/'S'ioo log LpjR log Lb logLx Area 

(Mpc) (Jy) (erg s —1 ) (erg s _1 ) (erg s _1 ) (arcmin 2 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



N253 


3.94 


n 


0.81 ± 0.15 


c 


1770.0 


1 


ISO 


0.67 


43.70 


42.80 


39.08 


y 


595.68 


N891 


9.08 




77 +o.i3 




126.0 


in 


ISO 


0.40 


43.69 


42.67 


39.78 


z 


143.14 


N1291 


8.6 


q 


1.07 ±0.15 


c 


10.1 


d 


IRAS 


0.17 


42.20 


42.74 


39.14 


f 


74.47 


N1482 


22.0 


a 






45.8 


a 


ISO 


0.77 


44.0 


42.42 


40.57 


z 


5.06 


N3184 


11.12 


r 


1.11 ± 0.22 


c 


29.0 


c 


IRAS 


0.31 


42.77 


42.33 


39.18 


ab 


42.24 


N3628 


7.91 


s 


0.82 ±0.20 : 


b 


106.0 


c 


IRAS 


0.49 


43.56 


42.50 


39.74 


z 


172.03 


N4214 


2.94 


o 


0.76 ± 0.20 


g 


29.0 


c 


IRAS 


0.62 


41.98 


40.89 


38.53 


g 


59.52 


N4244 


4.49 


1) 


0.46 ± 0.30 : 


b 


208.7 


c 


ISO 


0.40 


41.65 


41.35 


37.85 


z 


216.95 


N4449 


4.21 


1) 


0.51 ±0.20 : 


i 


73.0 


h 


IRAS 


0.49 


43.05 


41.48 


39.23 


i 


37.82 


N4945 


3.60 


p 


0.70 ±0.20 : 


b 


1415.5 


d 


IRAS 


0.42 


44.40 


42.45 


39.53 


z 


93.65 


N5253 


3.14 


t 


54+ - 21 


k 


29.4 


j 


IRAS 


1.04 


41.94 


40.81 


38.35 


k 


19.79 


N6503 


5.27 


p 


0.37 ±0.30 : 


b 


28.9 


c 


ISO 


0.35 


42.38 


41.66 


38.85 


z 


20.47 


M82 


3.89 


u 


0.50 ±0.08 


e 


1351.1 


d 


IRAS 


0.94 


44.27 


42.45 


40.63 


z 


98.87 


M83 


4.51 


X 


1.38 ± 0.28 


e 


638.6 


d 


IRAS 


0.41 


43.65 


42.59 






127.32 


IC5332 


4.63 


V 


1.30 ±0.31 


e 


5.1 


d 


IRAS 


0.16 


41.87 


40.97 






39.98 


Antennae 


13.8 


w 


0.53 ± 0.07 


c 


82.0 


c 


IRAS 


0.59 


44.30 


42.54 


40.43 


aa. 


22.05 



NOTES.- Column (2): distance in Mpc, see text in §2 for distance discussion for NGC891. Column (3): modulus of slope 
generated from the luminosity function. Data with a : has had the error estimated. Column (4): ISO/IRAS 100 \xm flux in 
Jansky. Column (5): Ratio of IRAS/ISO 60 to 100 um flux ratio. Column (6): Far IR luminosity calculated using EqsEJ 
and 171 Column (7): K-band luminosity calculated using equation 151 Column (8): the X-ray luminosity (logLx) estimated 
in the 0.3-2.0 keV energy range. Column (9): Deprojccted Galaxy area = 7ra 2 , where a is the semi-major axis of the D25 

isophotal ellipse. 
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